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r N ' The low frequency tail of the CMB spectrum, down along the radio range (~1 GHz), may 

'^ ' carry weak spectral distortions which are fingerprints of processes occurred during different 

,.i^ , epochs of the thermal history of the Universe, from z ~ 3 x 10^ to reionization. TRIS and 

Mh' ARCADE2 are the most recent experiments dedicated to the exploration of this chapter of 

CMB cosmology. The level of instrumental accuracy they reached in the determination of the 
5_( . absolute sky temperature is such that the removal of galactic and extra-galactic contamination 

"vi ' is the true bottleneck towards the recovery of the cosmological signal. This will be certainly 

C^ , the case also for future experiments in the radio domain. Here we present an update of a 

study originally done to recognize the contribution of unresolved extra-galactic radio sources 
to the sky brightness measured by TRIS. Despite the specific context which originated our 
analysis, this is a study of general interest, improved by the inclusion of all the source counts 
available up-to-date from 150 MHz to 8.4 GHz. 

m 

^ \ 1 Scientific case 

in 

^? , The frequency spectrum of the Cosmic Microwave Background (CMB) recovered by the FIRAS 

— I ' istrumentL!J on board the COBE satellite is an almost perfect planckian spectrum at the ther- 

modynamic temperature of 2.728 it 0.004 K, even if a more recent analysis and an update of this 
, , result can be found in a recent work by Fixsenl^. No signatures of spectral distortions have been 

r^ , found on the monopole scale. Nevertheless, there are important arguments'^ suggesting that 

C^ \ deviations from a blackbody curve originated in the pre-recombination Universe, may be present 

in the CMB spectrum. In particular, we know that if an energy injection in the photon-baryon 
fluid occurred in the redshift range ~ 10^ < z <~ 3 x 10^, corresponding to a temperature of the 
photon-baryon fluid ~ O.lkeV < T <~ IkeV, than the system relaxes towards a kinetic equi- 
librium, the photons having Bose-Einstein spectrum with a chemical potential /x 7^ 0. Several 
mechanisms have been proposed as possible sources of perturbation, among them decay of relic 
massive particles'^, dissipation of primordial magnetic fields'^, annihilation of relic particles'^ 
and dissipation of acoustic oscillations in the fluid'^. Was the spectrum a Bose-Einstein one, we 
would observe a dip in the plot Temperature vs Frequency of the CMB at frequencies <~1 GHz, 
given the present estimate of Hq (Hubble parameter) and Vti, (baryon density). The amplitude 
of this distortion, AT/T ~ ^{Qbh'^)^'^^'^ {h being Hq in units of 100 (km/s)/Mpc), is strongly 
model dependent. Extrapolating to lower frequencies the FIRAS upper limit on //, we expect a 
temperature dip smaller than few tens of mK. 

Another effect able to produce a monopole scale spectral distortion is free-free emission 
during reionization. This effect produces a signal AT/T oc YjjX'^, which is essentially the 



optical depth to bremsstrahlung, A being the wavelength and Yjf a distortion parameter as 
defined by Bartlett and Stebbinsl^. Despite new tighter limits on Yjj obtained by Gervasi et 
aLl^, we are still far from the possibility of testing reionization scenarios like those investigated 
by Weller et al.M 

Both effects are more relevant at decimetric wavelengths, where (1) calibrations are more 
difficult due to the size of antennas, (2)the brightness of the Galaxy overcomes the CMB and 
(3) the signal of Unresolved Extragalactic Radio Sources (UERS), the subject of this study, may 
introduce a temperature offset if not properly evaluated and subtracted. 

Two recent experiments searched for low frequency CMB spectral distortions: TRIS - and 
ARCADE2 . In the framework of TRIS, a new estimate of the brightness temperature of 
UERS has been proposed by Gervasi et al. L!^. Here we present an update of those results, 
exploiting all the most recent radio source counts from 150 MHz to 8.4 GHz and adding a new 
frequency to those used in the previous study. 

2 Unresolved Extragalactic Radio Sources contamination 

2.1 Aims and Methods of this study 

As stated in the previous section, the purpose of this work is to develop the earlier work by 
Gervasi et al.^^ enlarging the data-set described in table 1 of this reference. Here we included 
all the deepest radio source counts done in the last few years, especially those obtained at 
VLA and GMRT. In particular, we enriched our analysis taking into account new results at 150 
MHzEl, 610 MHz 1^ and 1400 MBz^ISl. We completed the 8400 MHz counts by adding the 
results by Henkel and Partridge^L^. 

We added also a new frequency in our analysis, namely the 325 MHz channel (Oort et al.^^, 
Owen et ahi^^). 

Typical experiments looking for spectral distortions in the CMB spectrum have beams > 7° 
(e.g. FIRAS), so that, on average, pointing in different directions, the radiometer will detect 
the same blend of AGNs, quasars and normal galaxies, with only a poissonian fluctuation in 
their number. Therefore, UERS are seen as an isotropic diffuse radiation. We calculate its 
temperature in two steps. First, we fit the differential source counts normalized to the euclidean 
counts, that is Q{S) = S'^'^dN/dS, S being the flux. Then we use the definition of brightness 
temperature, 

'^b,UERS = 7T7- / Q{S)S-'/^dS (1) 
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[ks Boltzmann constant), to calculate the integrated contribution of sources between two limit- 
ing values of flux. Here a problem is evident: given the fact that a survey will be complete at a 
flux limit Smin, how can we take into account the contribution of sources fainter than that? If 
we stop our integration within the data range, we get a lower limit of the temperature of UERS, 
introducing a bias in our analysis. A solution is to extrapolate our integral at fainter fluxes. In 
absence of a physical cut-off, to circumvent this problem we simply look for a functional form 
of Q(S) such that Tb^uERS remains finite when Smin — > 0, avoiding a kind of Olbers'paradox. 

2.2 Radio Source counts 

As originally proposed by Gervasi et al.'^we assume 

Q(S) = «,(5) + %(5) = -j^^,,l^^ + ^j^4^^ (2) 



where Ai,Bi,£i and /3j {i = 1,2) are parameters to be fitted. This analytical description, inspired 
by source evolutionary models proposed by Danese, Pranceschini and collaborators (see Danese 
et al.l^ and Franceschini et al.^^^, has the property of being integrable at faint fluxes, even 
if it is not good as log(S) polynomials to describe features within the experimental range. One 
of the basic assumptions underlying our approach is that the four parameters of Q{S) giving 
the slopes (ej and /3i) are frequency independent. Moreover, following again a suggestion of the 
model by Pranceschini et al., and analyzing the data at 600, 1400 and 5000 MHz (frequencies 
where both faint and strong fluxes are well sampled), we fix £i = £2, simplifying further the 
description of source counts. Now, the Qi and Q2 terms, following the notation of Gervasi et 
al.^^, are dominant respectively in the strong and faint flux regimes, so that two ratios can be 
calculated at the three frequencies with the widest flux coverage: r^ = A2/A1 and rs = B2/B1, 
at 600, 1400 and 5000 MHz. These ratios, r^ in particular, can be used to reconstruct the faint 
flux tails of counts at those frequencies where data are not deep enough. Since they show a 
weak frequency dependence (see tab. [T]for A2/A1), our choice has been to use the r^ and r^ 
calculated at the frequency closest to the one we want to reconstruct. 

Table 1: The measured A2/A1 ratios. 

I 600 MHz I 1400 MHz I 5000 MHz 



A2/A1 0.17 ±0.02 0.23 ±0.01 0.31 ± 0.03 



3 Results 

The main result we have obtained is a new estimate of the brightness temperature of the blend of 
UERS in the range 150 MHz - 8400 MHz. The calculated temperature values are well described 
by a single power law (eq. [3j) with spectral index -2.75 (see figH]). Our estimates are well in 
agreement with those found by other authors'^, and the main result is 

/ U \ -(2.75±0.02) 

T,^^«5(.) = (0.91 ± 0.02) (^^^^^^) K (3) 

for the contribution of unresolved sources to the overall temperature of the radio sky. 

This formula allows us to calculate the UERS temperature in correspondence of TRIS fre- 
quency channels, namely 600, 820 and 2500 MHz: we find, respectively, 950 ± 20mK, 408 ± 9 
mK and 20.6 ± 0.7 mK. 

4 Conclusions 

This update of our previous work allowed us a better reconstruction of differential source counts, 
especially at faint fluxes. In fact, unlike the first study, we could determine the ratio A2/A1 
finding a hint for weak frequency dependence. Improved source count fits lead directly to a more 
accurate estimate of the brightness temperature: this is very well described by a single power- 
law frequency spectrum. Then, using eq. [3l we could evaluate the UERS contribution at TRIS 
frequencies with uncertainties almost a factor 10 smaller than those commonly assumed before 
new estimates by Gervasi et al. "^. By virtue of that, the subtraction of this foreground is no 
longer a limiting factor in low frequency CMB experiments, the overall error budget being now 
dominated by the reconstruction of the absolute temperature scale and by the disentangling 
of galactic and cosmological signals. Finally, we conclude that the estimate given by eqlSj 
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Figure 1: Left panel: Temperature vs Frequency plot; calculated temperatures and power law fit (solid line); la 
error bars. Right panel: the same as in left panel, but with temperature rescaled by a factor oc i/-^''^ to appreciate 

the scatter of points around the single power law fit. 



exclusively based on source counts, doesn't agree with the extragalactic radio excess detected 
as a diffuse signal by the ARCADE2 collaboration^^. The discrepancy is around a factor 5 
at 600 MHz, and still the origin of the radio excess remains mysterious. If its origin is truly 
extragalactic, there must be a number of sources able to produce an intense integrated signal, 
but individually so weak to escape from detection with radio interferometers. The case is 
intriguing and again tells us how non trivial is searching for CMB spectral distortions, and how 
unpredictable are the outcomes. 
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